Entanglement between charge qubits induced by a common dissipative environment 
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We study entanglement generation between two charge qubits due to the strong coupling with a 
common bosonic environment (Ohmic bath). The coupling to the boson bath is a source of both 
quantum noise (leading to decoherence) and an indirect interaction between qubits. As a result, 
two effects compete as a function of the coupling strength with the bath: entanglement generation 
and charge localization induced by the bath. These two competing effects lead to a non-monotonic 
behavior of the concurrence as a function of the coupling strength with the bath. As an application, 
we present results for charge qubits based on double quantum dots. 

PACS numbers: 03.65.Ud, 03.67.Mn, 73. 21. La, 73.40.Gk, 73.63.Kv, 85.35.Be 
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I. INTRODUCTION 

Solid state nanostructures have become promising can- 
didates for quantum information processing 1 , with ba- 
sic operations like single-qubit manipulation and read- 
out having been demonstrated during the last few years. 
However, to go beyond single-qubit manipulations, and 
study effects such as entanglement generation and quan- 
tum gate operations, one needs some kind of interaction 
between the qubits. Although this interaction usually 
comes from a direct coupling between qubits (like the 
Coulomb interaction for charge qubits or exchange cou- 
pling for spin qubits), entanglement can be also gener- 
ated by coupling two qubits (which do not interact with 
each other) to a common third sy st en^iMi^iLS fl 1 1 1 1 . 
In most of these studies, the indirect interaction comes 
from the coupling to one or a few external degrees of 
freedom. Examples include the coupling to electromag- 
netic modes in a cavity (see, for example, Ref. 0, where 
the authors study entanglement of atoms within a single- 
mode cavity field) or to a harmonic oscillator represent- 
ing a mode in a thermal environmen t 10 ! 11 . Importantly, 
entanglement can be also induced when the environ- 
ment is made by an infinitely large number of degrees 
of freedom, namely a bath, as demonstrated by Braun 
in Ref. [Hj. This is an important case because entangle- 
ment is generated exclusively by incoherent means. In 
this context, different works have studied the coupling 
of two non-interacting qubits to fermioni o 13 i 14 ' 15 i 16 i 17 or 
bosonioii^OS^^ bathg 

Indirect qubit interactions have attracted attention 
recently because the information distribution among 
distant entangled particles is the base of quantum 



cryptography , quantum teleportation 
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dense cod o 24 i 25 ' 26 , different processes proposed for test- 
ing Bell inequalitie o 27 ' 28 ! 29 ' 30 and even certain steps 
within quantum computation algorithms. The pos- 
sibility of entangling two quantum systems which do 
not interact directly is therefore highly desirable, with 
various aspects of current interest like "entanglement 
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FIG. 1: Two independent double quantum dots coupled to 
a common phonon environment. Interdot hopping, t c , is al- 
lowed only in each double dot. The first double quantum dot 
is open to electron reservoirs, with probabilities for tunneling 
in and out given by Tl and Fr, respectively. 



swapping" 27,31 ,32^3,34 an d "entanglement transfer"^ 3 ^ 
In this paper we study entanglement generation be- 
tween two charge qubits due to the strong coupling with 
a common bosonic environment (Ohmic bath). For con- 
creteness, we focus on charge qubits based on double 
quantum dots (DQDs) but we point out, that our results 
can also be applied to Cooper Pair Boxes in a resistive 
environment. In a DQD, the electron charge degree of 
freedom is used to construct a qubi t 37 ' 38 i 39 ' 40 i 41 ' 42 , with 
logical states |0) and |1) corresponding to the localiza- 
tion of one excess electron on each one of the quantum 
dots (QD). One of the advantages of these charge qubits 
is their controllability through external voltages han- 
dling, as demonstrated in recent experiments^ 9 - where the 
charge has been coherently manipulated. We model the 
two DQD system as two independent two-level systems 
strongly coupled to the same Ohmic bath (two spin-boson 
models). In addition, we consider that one of the DQDs is 
coupled to electron reservoirs^ 3 - (Fig. []}, in order to allow 
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electronic transport. The coupling to electronic reser- 
voirs is treated using a Markovian approac h 42 ' 43 i 44 ' 45 , 
which is valid in the sequential tunneling limit and large 
bias voltages. The non-Markovian character of the strong 
coupling with the boson bath is, on the other hand, taken 
into account by using a polaron approac h 42 ! 43 ' 44 ! 46 ' 47 . As 
a result of this strong coupling, an indirect Ising-like in- 
teraction between qubits is induced by the bath. 

By combining both Markovian and non-Markovian ap- 
proximations, we derive a master equation for the re- 
duced density matrix of the system, including boson cor- 
relation functions in Laplace space. The resulting den- 
sity matrix is used to calculate the degree of entangle- 
ment (quantified by Wotters' concurrence 4 ^) as well as 
the probability for each one of the Bell states as a func- 
tion of the coupling strength by the bath. 

Our results complement previous work by Vorrath and 
Brandes in Ref. ,20j who studied a similar problem within 
a Markovian approach. We also mention some recent 
work s 11 ' 18 ' 19 in which related models are treated. 

The paper is organized as follows: in section II the 
model describing the DQD coupled to both, electronic 
reservoirs and the bosonic environment is discussed. We 
also present in Section II the general solution scheme 
for the density matrix equations. The coupling with the 
leads is treated by using a standard Born-Markov ap- 
proximation whereas the strong coupling with the bath 
is treated within a polaron approach. Section III shows 
the main results obtained, and finally we conclude in sec- 
tion IV. 



II. MODEL 

An array of two parallel DQDs in the strong Coulomb 
Blockade regime, and coupled to the same bosonic en- 
vironment, is considered. Interdot tunneling is allowed 
only in each double dot, defining an array of two charge 
qubits (Fig. [I}. The first DQD is weakly coupled to two 
electron reservoirs in such way that electronic transport 
through this double dot is possible (the excess charge in 
this DQD fluctuates between one and zero). The second 
DQD is closed and, therefore, has always one excess elec- 
tron. Note that such configuration is close to the one 
realized in very recent experiments^.The Hilbert space 
includes two-particles states |1) = \L1L2), |2) = Lii?2), 
3) = \R1L2) and |4) = I-R1-R2} (where Li(Ri) represents 
the charge localized in the upper (lower) QD of the i-th 
DQD), as well as one particle states |5) = IO1L2) and 
1 6) = IO1-R2} (where means no extra electron in the 
first DQD). The completeness of the system is therefore 

given by ^Li l fc )( fc l- 

The total Hamiltonian describing this system reads 

H = H sys + H res + Hb + Hsr + Hsb- (1) 

The free part of the Hamiltonian, i.e without couplings, 
contains three terms. The first term corresponds to the 



Hamiltonian of two independent DQDs, which in pseudo- 
spin language can be written as 



Hsys =53-A Ei (7j+t c <4, (2) 

i 

where Ae^ is the energy difference between quantum dots 
of each pair being Aei = el x — £Ri and A^2 = £l 2 ~ 
e r 2 (gj . is the on-site energy of the i-th QD of the pair 
j), (Tj is the j-th Pauli matrix acting on each DQD, t c 
is the electron tunneling amplitude which is considered 
identical for both DQDsi 3 -^. 

The Hamiltonian of the reservoirs, referred as "L" and 

Hres = ^{4 c L L C kX + ^k C k,R C k,R}^ ( 3 ) 
k 

where c\ ^ and Ck,p are fermion creation and annihila- 
tion operators in lead (3 with corresponding energy e^. 
Finally, the third term corresponds to the boson bath, 
which is described as a set of harmonic oscillators with 
frequency lu q : 

Hb = ^ Tiu q a\a q (4) 

a^ q (a q ) is the annihilation (creation) boson operato r 43 ' 50 . 

As we have mentioned already, we take into account 
the coupling of the system to both electronics reservoirs 
and a boson bath. The first coupling is given by 

H SR = { V k(4,L s LA + c-c) + Vk R (Ck,R s R,i + c - c -)}> 

fe,i£l,2 

(5) 

with Vl being the coupling with the lead (3. The 
Lindblad-type operators sla {sr,%) describe tunneling 
into (out of) the first DQD taking into account the 
two possible configurations Q^ ie i 2) m ^ ne secon d DQD, 
namely sl,i = |5)(1|, s ij2 = |6)(2|, s R:1 = |5)(3| and 
sr.2 — |6)(4|. Although we consider the coupling of only 
one of the DQDs to reservoirs, the generalization for both 
double dots is straightforward. 

The electron-boson interaction is described with a 
spin-boson Hamiltonian, where the bath "force" oper- 
ator = J2 q lq( a l + a q) couples linearly to each qubit' s 
g^ . 50 ' 51 (here we consider that such interaction is identical 
for both DQDs, 52 7* = 7,): 

i7sB = EE^K+ a 9)> ( 6 ) 

q i 

Note that this coupling which is longitudinal in the lo- 
cal basis of each qubit, contains both longitudinal and 
traversal components in the basis that diagonalizes the 
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qubit Hamiltonian 5 -. The bath effects can be encapsu- 
lated in the spectral density J(u) = ^2 q l 2 &{w — oJ q ). 
In the following, we use a generic ohmic bath: J(u>) — 
2awe" u ^ c , where ui c is a cutoff frequency and a is a 
dimensionless parameter which reflects the dissipation 
strengt h 50 i 51 i 54 . As we shall see in the next subsection, 
the coupling of both qubits to the same quantum heat 
bath leads to decoherence and to an effective interaction 
between qubits. 

A. Polaron Transformation 

Due to the strong coupling of the qubits with the bo- 
son bath, a proper description of the system must take 
into account non-Markovian effects. Among the differ- 
ent approaches available to deal with thi a 55 i 56 ' 57 i 58 , we 
use a "polaron transformation"—, which for an arbitrary 
operator O is given by: 

0= e s Oe~ s , 

s = EEr?K-4 (7) 

q i 1 

This approach, which is well-known for treating problems 
in which bosonic modes couple to localized electronic 
states, has been successfully used for studying single 59,60 
and double quantum dot o 42 i 43 ' 44 > 46 strongly coupled to a 
bath of phonons. 

Applying the transformation in Eq. ([7]) to the rele- 
vant operators in our model, we obtain the transformed 
operators: 





< 


K = 




a q = 


H % 


SL,t = 




SR,i = 





where a\ and a l _ are the ladder spin operators on each 
DQD. The operators X and X^ are polaronic phases^ 
given by 

X = e A , (9) 

with 

A = ^ (4 ~ a i) ■ 
i 

(10) 



By substituting these transformed operators into the 
equations ((2|) to (|6|) , we obtain the effective Hamiltonian: 

H= H + H T + H SR (11) 
H = ]T \ Ae t al « + H B + H res (12) 



The effect of the canonical transformation is threefold: 

(i) The electron-boson interaction H$ b has been trans- 
formed away. 

(ii) The state of the bosonic system is strongly modi- 
fied every time an electron tunnels between dots (boson 
" shake-up" As a result, the interdot tunneling ampli- 
tude (Eq. I13p becomes renormalized with environment- 
dependent phases through the operators X = e A . These 
time-dependent exponential phases, which appear as a 
result of the non-perturbative treatment of the electron- 
boson interaction, lead to non-trivial effects. In particu- 
lar, this implies that non-Markovian effects become rel- 
evant and need to be considered in the dynamics of the 
reduced density matrix. Note also that, in principle, this 
renormalization of tunneling has to be taken into account 
also in H $r through the operators ~s~L,i = SL,i e ~ A ^ 2 an d 
sra = sji.ie A l 2 . However, this is no longer true in the 
limit of large bias voltages, where the coupling to the 
reservoirs becomes Markovian (see the next subsection). 

(iii) The transformed Hamiltonian contains an effec- 
tive interaction between qubits H e ff = —j^2ijO~lal 
due to the coupling with the common bath; this inter- 
action has an Ising form and depends on the parameter 
k = ^2 q ~fq/oj q (which for Ohmic dissipation used here 
reads n = 2auj c ) favoring states with the same charge 
distribution in both DQDs. 



B. Master equation 

We define the total density operator of the open sys- 
tem as x(t) — e~ lHt x(0)e lHt which, after the transforma- 
tion in Eq. ([7]), can be written in the interaction picture 

as 42 i 44 ~ = e iH t^ e -iH t^ with ^ = e -iHt^^ e iHt 

By taking the partial trace over the reservoir degrees 
of freedom, the reduced density matrix (RDM) of the 
two DQDs plus the boson bath is obtained as p(t) = 
Tr res x(t). Applying the second order Born approxima- 
tion we obtain the equation of motion for p(t) as: 
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Jt p(t)= -i[H T (t),p(t)} 



E 



k,ie~L,2,jEL,R 



E 



k,i£l,2,j£L,R 



E 



k,i£l,2,j£L,R 



df \vi\ 2 f j (4)e^ t - t 'Hs j MAi(t')p(t') - AS'W^m 



dt' |V^| 2 [1 - / J (4)]^ teJfc( ^*' ) {^, J Ci)^, i (i / )^') - ^^(iOp^)^^*)} 



df \vi\ 2 fHiy^'Hms^Ut) - AiVw^itf)} 



E 



k,i£l,2,j£L,R 



r 



(14) 



where = Tr res {i?oc|, jCk.j} are the Fermi dis- rates in and out of the DQD. Working in an "infinite 

tributions of each contact (i? is the density ma- bias regime" between the reservoirs (such that f L -> 1 

trix of the electron reservoirs, considered in ther- and f R — > 0) and assuming a constant density of states 

mal equilibrium) 42,44,46 Eq ^ can be sim- in the reservoirs, the coupling with the leads becomes 

plified by rewriting the sums over k as integrals Markovian: ~J2 k \Vfc\ 2 f L (ei)e le >> lyt ~ t ) =Ti,5(t — t') and 

E k \Vi\ 2 P(4)^- t,] = IZo #r j -( e )^( e )e«(*-*'), £ fe l^| 2 [l " / fl (ef )]e^(*-*') = T R S(t - t<) and, there- 

where T (e) = 2nY,k \ V k\ 2s ( e ~ 4) arc thc tunneling fore > E q- d reads 

I 



d ~, x 



-i[H T (t),p(t)) 

~ E &R,iV)*n,iV)p(f) - 2^(W)4,(0 + P(*')Sj ll i(* / )S'fl,i(*')} 



(15) 



r 



As we mentioned already, the fact that the coupling with 
the reservoirs becomes Markovian in this limit implies, 
in particular, that the renormalization of tunneling due 
to the bosonic bath becomes ineffective (for example, 

^,i(O4 ) i(*0 = ^,i(O4,i (*'))■ 

Invariance under unitary operations implies that 
the expected value of any dot operator can be 
written as (0(f)) = Tr dot {Tr B {p(t)}Tr B {0(t)}} = 
Tr dotyB {p(t)0(t)} , where Tr B is the trace over the bath 
states. In particular, the expected value of the projec- 



tor operators over the system states Y nm = \ri)(m\, can 
be written as (Y nm ) = Tr dot {p s Y nm } = (m\p s \n) = 
Tr dot {p s Y nm }, where we have defined the RDM of the 
DQDs array (system) as p s = Tr B p. It is therefore pos- 
sible to obtain matrix elements of the reduced density op- 
erator by just calculating the expectation value for the 
suitable Y nm (t) operators directly from the master equa- 
tion P3| . Using the notation p^ nn = (m\p s \n), we obtain 
the following set of exact equations: 



J 



Pnm(t) = P„ m (0) - * / Tr dot , B {Kt')[Ynm(t), H T (t')}}dt' 

Jo 

/ Tr dot M(S L (t')si(t')p(t') - 2s{(t')p(t')sL(t') + p(t')Mt')sl(t'))Yrnn(t)}dt' 
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/ Tr dot M(sUt')s R (t')p(t') 
Jo 



The full expression of the density matrix is too large 
to give it here (its total dimension is 6 x 6) and we show 



Pis(*) = 

pUt) = 



with s n \n) = [J2i - 3 K E 4li a >i) \ n )- 

Note that Eqs. (|17[) are not closed. They con- 
tain expectation values involving products of dot and 
boson operators, as for example (Y 33 (t')X}Xt>) = 
Trdot,B{p(t')Y 33 (t')X}Xt>}, which need to be decou- 
pled. If one is not interested in the system back ac- 
tion on the bath, the latter can be assumed to re- 
main at thermal equilibrium at all times 42. Therefore, 
the reduced density operator can be approximated as 
p(t') ~ Pb(0) S5 TrBp(t'). By using this approxima- 
tion, we can decouple higher order correlation functions 
as (Y 33 (t')x}x t ,) « (Y 33 (t'))(x}x t ,), etc. This decou- 
pling corresponds to the so-called Non Interacting Blip 
approximation in the spin-boson problem^. 

For an equilibrium boson bath one can write the corre- 
lation functions asi 2 . C(t-t') = {X t x\,) = e -*(*-*') with 
$(r) = r °° ^{(1 - cos(wr))coth(/3w/2) + isin(wr)}. 
Note that in our problem we also need C^it — t') — 
(XtX t Xj,xl,) = e _2 *( t_t ) 5 which appear from coher- 
ences involving interdot processes like |1) <-> |4) (namely, 
\L\,L2) <-> |i?i,i?2), see Eq. (fT7|) ). In principle, terms 
involving half phases e~2*( t-t ) also appear in tunneling 

processes to the reservoirs (like, for example, \R\, L2) -? 
|0i, L2)) but, again, they do not contribute in the Marko- 
vian limit. 

The resulting set of coupled equations can be given in 
matrix form as 

Ps(t) = Ps (0) + f\M(t - t') Ps (t') + T}dt' (18) 
Jo 

where ps is a vector containing the different matrix el- 
ements of the reduced density operator, the vector T 
contains the terms related with the coupling of the first 
DQD to the reservoirs, and M(i — £') is a non-Markovian 
time dependent kernel which contains the bath corre- 
lation functions^ 3 . Equation (Ti~8"|) can be solved in the 
Laplace spac e 42,43 as 

p s (z) = [z- zM(z)]- 1 ( jOs (0) + T/z). (19) 



- 2s fl (i>(0^') + p(fp R (t')s R (f))Y mn (t)}df. (16) 

I 

just two examples for pf^t) and pf 3 (t) elements: 



I 

The kernel M(z) contains the Laplace trans- 
form of the bath correlation functions = 
J °° e~ ZT e^ tST C^ (t)(1t evaluated at different en- 
ergies corresponding to the involved transitiom 42 i 43 ' 44 

III. ENTANGLEMENT 

The full time-dependent density matrix can be ob- 
tained by solving algebraically Eq. (fiU)) and performing 
an inverse Laplace transformation, which is a formidable 
task. Fortunately, the entanglement generated by the 
bath is finite at long times, namely in the stationary 
state, as we will show. The stationary solution of Eqs. 
(1181) . poc, is obtained by extracting the 1/z coefficient in 
a Laurent series of ps(z) for z — > i 42 ' 43 i 44 For entangle- 
ment quantification we use Wootters' concurrence 4 ^ for 
a general state of two qubits, 

C = max{0 7 y/% - ^f\~ 2 - V^3 - \J%) (20) 

where the As are the eigenvalues in decreasing order of 
the non hermitian matrix ps(o~ y <S> o~y)p*s(°~y ® &y)- The 
concurrence ranges from C = for non-entangled states 
to C = 1 for the maximum degree of entanglement. 
That maximum entanglement is showed by the Bell 
states 48 . In the basis of triplet and singlet states. So) = 
^(liiifc) ~ \R1L2)), \T ) = ^(|Liik) + \RiL 2 )), 
\T + ) = IL1L2} and |T_) = \RiR 2 ), the Bell states read: 
|*+) = \T ), |*-) = |5 ), \<i> + ) - ^(|T_) + \T+)) and 

l^> = 72(l T ->-l T +»- 

Importantly, the stationary density matrix in our prob- 
lem corresponds to a transport situation and, there- 
fore, a proper generalization of concurrence to non- 
equilibrium is needed. Following Ref. UtI we quantify 
non-equilibrium entanglement via the concurrence C of 
the stationary state Ppca, where P is the projection 
onto doubly occupied states including proper normaliza- 
tion. The projection P corresponds to taking the limit 



-i{t c e- i ^-^ t - t '^(Y 33 (t')xlx t ,) - {yu(t)XtrX}) + (Y 32 (t'^xlx t ,) - (^(i')*^)) 
-i{t c e- i ^-^ t - t '\(Y^)Xt'X}x}x t ,) + {Y 2 \{t>)X t ,X t ,x\x t ,) 

+{Y^)X v x\x\X v ) - {Y 3 \{t')X t ,X v xlxl)), (17) 
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FIG. 2: (Color online) a) Concurrence as a function of the 
strength of dissipation a. b) Population of triplet and singlet 
states. Parameters: Aei = Ae2 = 0, t c = 3.5, Tl = 10, 
and LOc = 500 (in units of Yr — lfieV). These parameters 
correspond to typical experimental values in AlGaAs-GaAs 
lateral DQDs^^^. 



Tl — > oo where both qubits are always occupied with 
one single electron. For concreteness, we focus on the 
zero-temperature case. 

The concurrence as a function of the coupling a al- 
ways shows the same qualitative behavior: for very small 
a there is a switching behavior, indicating that below 
a minimum interaction strength k the concurrence van- 
ishes, cf. Fig. [5{a) for identical QDs (Aei = Ae 2 = 0). 
As a increases, two effects compete: entanglement gen- 
eration and localization induced by the bath. At small 
a, the two delocalized states \Sq) and \Tq) have a finite 
weight which depends in a nontrivial way on the ratio 
On the other hand, for strong coupling the bath com- 
pletely freezes the charges on the left dots and the triplet 
7+) = \Li,L2) becomes fully occupied, cf. Fig. EJb). 
These two competing effects lead to the non-monotonic 
behavior of the concurrence vs. a, with an optimal value 
at which the concurrence presents a maximum. 

The population of each Bell state is shown in Fig. [3] 
The system does not originate a preferred Bell state and 
therefore both maxima in the concurrence contain con- 
tributions from all states. The first concurrence peak is 
formed by a combination of the four Bell states with a 
symmetric contribution of | V P + ) and \(f> + ), whereas on the 
second peak \<f>~) probability is slightly dominant. Elec- 



0.5 



0.4 



S 0.3 
!5 

CO 

■g 0.2 



0.1 - 



0.0 













' . y ' 


p 


- 


/ 

.' s ' . •■ s 
< t' . . V 


W + 

P 

w- 

p 


- 


i ~f — ^\ N 

\ .' ^ 


p 













0.00 0.05 0.10 0.15 0.20 0.25 0.30 
a 

FIG. 3: (Color online) Population of the Bell states as a func- 
tion of q. Same parameters as Fig. [2] 
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FIG. 4: (Color online) Color map of concurrence vs. interdot 
tunneling t c and a. The rest of parameters are the same as 
in Fig. H 



trons localization in "parallel" charge states is reflected 
in the large probability for both \cf> + ) and \cf>~) states for 
a > 0.2. 

The concurrence as a function of both t c and a is shown 
in Fig. fll For 2t c < Tr, the dephasing induced by the 
leads suppresses interdot coherence and the contribution 
of the delocalized states \So) and \T ) is negligible. Thus, 
the concurrence is almost zero for all a. For 2t c > Tr, 
entanglement is hnite in a region a m i n < a < a ma x', 
both, a m i n and a max increase with t c . For 2t c » T/?, 
the system present a maximum in the concurrence at 
a w 0.15 with values C w 0.3. 

The effect of on concurrence is shown in Fig. 
[5] Here, we also find the switching behavior described 
above: starting from a = 0, the state of the system is 
strongly mixed for small T^. Therefore, C = below a 
minimal value a m i n . This threshold value decreases as 
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T/j increases. At fixed a, the entanglement decreases as 
one increases Tr. For very large r^, the pure localized 
triplet |T + ) = \Li,Li2) is reached and thus the entangle- 
ment is zero. This effect, which is a transport version of 
the Quantum Zeno effect, is similar to the one occurring 
in capacitively coupled charge qubits open to reservoirs ^ 

A finite detuning Ae^ > (Ae, < 0) localizes the 
charge on the lower (upper) QD of each pair and, there- 
fore, the entanglement should depend on whether Aei = 
Ae 2 > or Aei = — Ae 2 > 0. The concurrence of the 
latter case is very similar to the one for Aei = Ae 2 = 0, 
and therefore the population of singlet and triplet states 
show also the same kind of behavior, Fig. E£a). On the 
contrary, the concurrence for Aei = Ae 2 > is different 
with a narrow resonance at small a, cf. Fig. [5] This 
resonance corresponds to a maximum in the population 
of the triplet T , cf. Fig. [TJb), followed by a fast decay 
of both To and So and an enhanced population of T + 
(and, hence, zero concurrence). The overall qualitative 
behavior is in agreement with Ref. where the current 
through two DQDs coupled to the same phonon bath is 
analyzed in the Markovian limit. The indirect interac- 
tion due to bath leads to an enhancement of the inelastic 
current at Aei = Ae 2 > 0, and a maximum population 
of the triplet To, which is a transport version of the Dickc 
effect. In close analogy with the Dicke effect in quantum 
optics, this superrradiance should turn into subradiance 
as the probability of finding the system in the singlet Sq, 
rather than in the triplet To, increases^ We do not find, 
however, the subradiance counterpart in our analysis of 
concurrence. 
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FIG. 7: (Color online) a) Population of triplet and singlet 
states for Aei = Ae2 = 10. b) The same for Aei = Ae2 = 
-10. 



IV. CONCLUSION 

We have shown that the strong coupling of two in- 
dependent qubits with a common bosonic bath at zero 
temperature originates entanglement between the qubits 
in the stationary limit. We also identify that two ef- 
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fects compete as a function of the coupling strength with 
the bath: entanglement generation and charge localiza- 
tion induced by the bath. These effects lead to a non- 
monotonic behavior of the concurrence as a function of 
the coupling strength with the phonons. 

In addition, the concurrence strongly depends on tun- 
neling and energy difference on each DQD as well as on 
the coupling with external leads, parameters which can 
be controlled experimentally. 

Due to the small concurrence values obtained here 
(C < 0.5), and the fact that no preferred Bell state is 
formed, the use of this setup may not be an optimal 
choice for entanglement studies in the solid state realm. 
Note, however, that this system is the minimal imple- 
mentation of a fully tunable two qubit system coupled 



to a common bath. From this point of view, this real- 
ization is an attractive benchmark in which to study the 
interplay between quantum coherence, entanglement and 
decoherence. 
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